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Abstract 

The zero-range potentials of the radial Schrodinger equation are investigated from a point of Darboux transformations 
scheme. The dressing procedure is realized as a sequence of Darboux transformations in a way similar to that used to 
obtain the generalized zero-range potentials of Huang-Derevianko by specific choice of a family of parameters. In the 
present approach we stay within the framework of conventional zero-range potential method whilst the potential parameter 
(scattering length) is modified taken into account spectral molecular properties. This allows to introduce molecular zero-range 
potential once the corresponding discrete spectrum is known. The results are illustrated on example of flat cyclic molecular 
structures, with particular focus on a benzene molecule, which bounded states energies are first found using atomic zero-range 
potentials, compared with the Huckel method, and then used to introduce single zero-range potential describing the entire 
molecule. Reasonable scattering behavior for newly introduced potential gives a possibility to tackle many-molecule problems 
representing molecules as appropriate single zero-range potentials. 

1 Introduction 

Zero-range potential (ZRP) method is proved to be an efficient tool in describing quantum structures with arbitrary 
geometry at sufficiently low energies such that detailed structure and interactions of atoms (replaced with point- 
centers) between themselves can be neglected [T31 HJ [21 [SI [12] • Simplicity of the approach makes it very attractive 
and over the years it has been further developed into generalized zero-range potential (gZRP) method that takes 
into account effects of higher order partial waves in describing scattering properties [HI [TE[ [TU[ [H [TSl [IHl [IS]- 
However, the scope of application of the regular ZRP method can be extended to describe not only small, typically 
atomic, structures but larger molecular systems. 
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In the present work, after giving brief but self-consistent review of Darboux transformation (DT) and gZRP 
ideas, we enhance standard ZRP approach by means of appUcation of a dressing procedure [TTl [151 [HI ■ We notice 
and make use the fact that a certain class of DTs of spherically symmetric ZRPs yields also a potential of zero 
range but, in general, with different effective characteristic (s-wave scattering length of the original ZRP is altered). 
Transformation parameters of a single atomic potential can be chosen based on solution of conventional bounded 
state problem solved for a system of ZRPs (molecule) such that resulting solution mimics scattering behavior for 
the molecule as a whole at low energies. 

Described ZRP treatment is demonstrated on example of a generic cyclic structure with focus on a particular 
case of benzene molecule. 

2 Method of zero-range potentials 

The idea of the ZRP method is to replace real potential with a point potential placed in the origin that mimics 
physical scattering behavior in the far zone. 

Therefore the equation in question is the radial Schrodinger equation for the free space (r > 0) 



where = 2^ is wave number, /i, E are mass and energy of a particle (electron), respectively, h is the Planck's 
constant, -0 — ip^r) is the radial part of the wave function. 

It is well-known that the general solution can be formed as a linear combination of spherical Bessel and Neumann 
functions, jiikr) and yi(kr), respectively, 




(1) 



ip{r) = Ci {ji{kr) - t&nrii ■ yi{kr)) , 



(2) 



or as a combination of spherical Hankel functions 




(3) 



with si — exp {2irii) being a scattering matrix. 



Taking into account the following asymptotes at /cr — >■ [T] 



ji{kr) 



{kr) 



(4) 



{2l + l)\\ 



yi{kr) 



(2Z- 1)!! 
(fcr)'+i 



(5) 
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written with notion of the odd factorial {21 + 1)!! = (2Z + 1) • (2Z — 1) • . . . • 3 • 1, (—1)!! = 1, we obtain asymptotic 
behavior of finite-energy solution at the origin 



N ^ f (krY (2;-l)!!\ ^ (2Z-1)!! 



(2/ + 1) 



From here, the constant Ci can be expressed as 



(fcr) 



{kry 



' {21 + 1)1 dr^i+^ ^ ^' 



r=0 



Elimination of Ci yields 



1 



(2Z + 1)! 



_o (2;-M)!! (2Z-1)!! tanT?,^ 



This is to say that imposing the boundary condition 



1 d 



21+1 



2H\ 



(2^-1)!! 



1/a 



2Z-I-1 



with 



tan?7; 



(6) 



(7) 



(8) 



(9) 



(10) 



results in the desired physical behavior of solution at far distance from the point-center. 

Introduced quantity a/ is termed as the partial wave scattering length and is independent of particular energy 
value at low energies. To justify this fact, it is enough to consider zero-energy solution to the equation ([T]) 



1 



^p{r) = Air' +A2 — . 



By matching this with ([2]), we conclude 



tan 77; w 



(2Z-|-1)!!(2;-1)!! 



that demonstrates energy-independence of scattering length parameters a;. 

In most of the cases, at low energies, it is enough to consider ZRP as a spherical point, thus neglecting higher 
partial waves. Therefore, solution to the ([T]) should be a subject to the simple boundary condition 



dlog (r-f/') 
dr 



(11) 



r=0 
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where /? = 1 /oq is the s-wave inverse scattering length. 

3 Darboux transformation and its application 

Consider one-dimensional Sturm-Liouville equation 



= AV', (12) 

where 

Now we apply such algebraic-differential transformation (referred as Darboux transformation [8]) 

D = A-.(.) (14) 
that the transformed equation preserve the same form 



LW^w=A^w^ (15) 



where ij^-^^ = Dij)^ L^^l = — — r -1- v^^^x) and the transformed potential is 



u^^\x)=u{x)-2a'{x). (16) 

It is straightforward to check that 

rT(x^ — 

$1 

with $1 being a particular solution to the original equation, i.e. it satisfies 



$1 



+ u(x)$i = Ai$i 

for some Ai. 

It can be shown that in case of multiple application of the Darboux transformation, the expression (|17p is 
generalized into 

a(x) = -^[l0giy($i, ...,$Ar)] 

ax 

resulting in the following transformation of the potential 



3 Darboux transformation and its application 



5 



u^^\x) = uix) - 2— [logM^($i, . . . , $jv)] . (18) 
Corresponding transformation of the solution is given by the Crum's formula [71 [20] 

where denote determinants of the Wronskian matrices. 

Having considered one-dimensional Cartesian case, we move towards the radial Schrodinger equation ([ij and 
notice that it can be brought to the form eligible for direct application of obtained transformation formulas. 
Namely, performing substitution ip — x/i^^ one readily obtains 

-x"+^-^^X = k\. (20) 

That is to say, that we can apply Darboux transformation to the equation ([IJ meaning that all original wave 
functions ijj should be multiplied by r whereas the potential term 

= ^ (21) 

remains unchanged. 

Thus, the formulas (fT8|) . ((23|) in spherical coordinates should be replaced with 

u[^l(r) =M(r)-2^[logiy(r$i, r$w)] (22) 

and 

^ ^[^] ^ W^(r<&i,...,r$^.,rV.)^ 
^ ^ rVK(r$i, r$Ar) ^ ^ 

Following an idea of the previous works [llj, we first show how gZRP can be induced by application of DT. 
We start by choosing a spherical Bessel function as the seed solution 

Mr)=jiikr) (24) 

and apply A'^-th order Darboux transformation by taking spherical Hankel functions with specific parameters Hm 
as prop functions 

<^>mir)^h\'\~iK^r), m^l,...,N. (25) 
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We can employ Crum's formula (|19p and consider asymptotic behavior of spherical functions at r — ;> cxd 



ji{kr) 



sin {kr — Itt/2) 
kr ' 



(26) 



Vi{kr) 



cos {kr — Itt/2) 
kr ' 



(27) 



i+iexp(-ifcr) 



('21 

h\ '{kr) = ji{kr) - iyi{kr) w i 



kr 



Then the Wronskians turn into Vandcrmond determinants, hence, 



(28) 
(29) 



i>f'\r)^C 



exp {ikr) A (ki 



kr 



A (ki, . . . , Km) 



, kat, ik) .1 exp {—ikr) A (ki 



A . . . , hn) 



(30) 



where C represents here and later on a generic constant without prescribing it any specific value such that it can 
absorb constant multipliers where their meaning is not important without changing notation. 
We note that for I = this expression is not asymptotic, but exact. 

A Vandermond determinant in (j30p can be computed by noticing that k = —inm (for to = 1, ... , N) are the 
roots of polynomial with respect to k equation that is obvious due to the fact that replacement ik yields zero 

determinant due to linear dependencies of the rows, thereby allowing the following factorization 



A (ki, . . . , UN-, ik) 



Kl 
K2 



"-1 

K2 



1 Kjv K 



N 



ik {ikY 



N 



N 



C Y[ i'^m - ik) . 



Denoting 



si 



N 

n 



N 

n 



{k + iKm) 



{Kjyi + ik) {k iKm) 

we recognize in the expression (I30p the asymptotes of spherical Hankel functions, hence 



(31) 



i/,|^l(r) = C sihX\kr) - hY>{kr) 



(32) 
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As one can notice, the expression ([5^ coincides with ([3]) if 



exp(2z^,)=nr'^^' (33) 

[Km + IK) 



or, taking into account pil| . 



tanry, = -af+^0+^ = _.^ Tll=ii'^rn ik) H^^i Kn + ^fc) ^ ^^^^ 



Thus we conclude that for the direct correspondence to the gZRP behavior, the number of consequent transfor- 
mations = 2^ + 1 should be taken. 

Given some quantity a = \a\ e"^", the parameters ki, . . . , k,21+i needs to be chosen in a way that 

21 + 1 

Jl {Km + ik) = (ikf^^ - a. (35) 
This is equivalent to the ik = —Km {m = 1, . . . , 2Z + 1) being the roots of the equation 

(ik)^^^^ = a. 



that is to say —Km = or 



- Km = \a\ exp 



In the similar manner we obtain 

21 + 1 



l/(2/+l) / (0Q+27rTO)\ 



* 27Tl j' "^ = 1. •••>2/ + l. (36) 



Jl («:„, - ifc) = {-ikf+^ - a. (37) 

m— 1 

Substitution of ([25]), (I2Z1) into resuhs in 

a = (-l)'+Var^ 

Therefore, providing a; is a real number, p6p yields 

/ ; 4_ 2to + 1 \ 

Km = -l/fli ■ exp ( ITT — 2l-irl — )' TO = 1, ■ • ■ , 2/ + 1. (38) 

Now alternatively to gZRP instead of taking into account impact of higher harmonics, we consider the trans- 
formation of spherically symmetric solution (i.e. / ~ 0) with the parameters Km to be chosen according to spectral 
data of the entire molecule obtained, for example, from the discrete spectrum of molecule modeled by regular ZRP 
method. 
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We get back to (|34p and since 

N N N N N N N N 

Y[ (k™ + ifc) = Y\_ '^rn + X! n + X] 5Z n + • • • + ^ Kn + (j^)^ , 

m—1 m— 1 n— 1 m— 1 — 1 m— 1 n— 1 

we continue the last equality as 

«fc 11™=! '^m + . . . + {ik)^ 

^21+1^21+1 ^ ^ ^ . (39) 

rim^l Km + • • • + [ik) X]n=l '^n 

Limiting ourselves with I = case, we consider low-energy limit neglecting terms of higher order than one with 
respect to k to obtain 

A^n=l llm=l 

ao = ^ • (40) 

Alternatively, recalling notation /3 = I/oq, the last expression can be rewritten as 

(41) 



\m=l / 



On the other hand, we can readily observe that Darboux transformation of the seed solution with the 
prop functions (j25p results in the solution corresponding to the ZRP for arbitrary choice of N and parameters Km, 
m—1, . . . , N. Indeed, for / = 

N 

\ogW{r^i, r$Ar) = K,nr + C, 



thus, according to (|22|) . 

u[^l(r) = 0, r > 0. 

Freedom of choice of transformation parameters gives a possibility to induce desired poles of scattering matrix 
(|3ip and, thereby, perform transition from atomic to molecular ZRP by choosing m—1, . . . , M such that new 
potential allows M bounded states obtained by solving conventional formulation of the discrete spectrum problem 
for the set of M standard ZRPs making up the molecule. Molecular ZRP can be qualitatively characterized by 
effective scattering length computed according to PU)) and thus used in such simplified scattering calculations in 
complicated problems involving, for example, a chain of molecules. 

Below we demonstrate the presented approach by applying it to a general cyclic molecule with eventual focus 
on a benzene molecule. 
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4 Case study: a cyclic molecule 

In order to describe a cyclic structure, we consider a system (molecule) of M point-centers, typically atoms, modeled 
by ZRPs which positions are kept fixed in the plane in space in a cyclic order. 

Spatial coordinates are chosen such that all atoms lie in the xy-plane with the origin coinciding with the center 
of the molecule. Atoms are numerated counterclockwise and the x-axis is chosen along the direction to the first 
atom. Described geometry is illustrated on the Fig. [T] 

Z 

A 




X 

Fig. 1: Illustration of the geometry by an example of 3 atoms. 



4.1 Bounded state problem 

According to the idea of ZRP method, we write solution to the bounded state problem as 



M 



(42) 



subject to the conditions (fTTj) 



'log( 



r-Ri 



z = ...,M, 



\r-Ri\=0 



where 



9 (r, n) 



(43) 



(44) 



Let us introduce the following notation. We will denote a distance between nearest neighboring atoms as Ai?i , 
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between every second atom from each given one as Ai?2, every third as AR3, and so on (see the Fig. [5]). 




Fig. 2: Illustration of the distance notation. 



From simple geometrical considerations it follows that 



AR, = 2i?o sin ( ^ • I ) = 2Ro sin (0j+i/2) , 



(45) 



27r 



■j — — [j — 1) is azimuthal angular coordinate of a j-th atom, Rq is radius of the molecule. We note that 



where 

the expression is formally valid for j — 1, . . . , M, however we note that there is only [Af/2] different values since 
Ai?[^//2]-j — ARl^]^J/2]+j for j = 1, . . . , [M/2] — 1 where [•]. here and later on, marks integer part of an argument 
(i.e. the floor function). 

We develop conditions (|43|) to show existence of the derivatives 



9 log (^f~Ri ■ ilj{f)j 



d r- ft 



\r-Ri\=0 



[(. . . + C,_i • g (Ai?i, k) - C, • K + C,+i • g (Ai?i, n) + ...)- 



r - R, 



d {i}{r) -g(^ 


r — 


R. 


,.)) 


d 


r - 


-R^ 





\f-Ri\=0_ 

To show that the last term in the square bracketed expression is well-defined and equal to zero due to the 
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presence of the multiplier 



r-Ri 



we first note that ip{r) — g 



Rr 



Rj 



for J = 1, ... , N, j 7^ i, which does not lead to a singularity when computed at 

r-R 



chain rule, it remains to demonstrate that the derivatives 
in a straightforward manner as it follows 



r-Ri 



lij is a function of the arguments 

f—Ri =0. Thus, due to the 
f — Ri =0. This can be done 



do not blow up at 



d 


r — 


Rj 


d 


f— 


Ri + ARij 


d 


r — 


Ri 




d 


r-Ri 











f-R^ 



' + 2- (Ai?y , f-Ri^ 



ARij 



d 



f-Ri 



r- R, 



f-R^ 



(I (A/?,j. r- Ri 



d 



f-R^ 



where ARij = R^ — Rj. 



Taking into account that 

d(^ARij, r-R^ 



r-Ri 



lim 



(ARij, tj 



\P-Ri\=0 



AR, 



cos a 



where cos ay = l^^^j fj\ , we conclude 



d 



d 



r-R, 



r-Ri 



= cos a 



I] 1 



thus, we have shown the required boundedness 



-1 < 



d 


r — 


Rj 


d 


r — 


Ri 



< 1. 



\f-Ri\=0 



Therefore, eventually. 



o>log ( 


r- B, 


■ ' '(F)) 


d 


r-R^ 





— (... + Ci_i • g {ARi, K)-Ci-K + d+i ■ g (Ai?i, «) + ...) = -/3 

^i 



\r-Ri\=0 



for i = 1, . . . , M. 

This gives the set of homogeneous linear algebraic equations allowing to determine Ci, . . . , Cm while its com- 
patibility condition yields spectral equation with respect to k. 

Alternatively, the set of constants can be found from general symmetrical considerations and spectral condition 
will follow from their direct substitution into the last equation taken, for instance, at i = 1. Thus the spectral 
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condition reads 

M-l 

(3-K+J2Cj+i/Ci-g{ARj,K) = 0. (46) 

To find all constants Cj employing symmetry of the problem, we introduce cyclic permutation operator T 
which action on the wave function (|^^ results in the replacement Ri — > R2, R2 — > -R3, Rm Ri- Since 
T^'^ilj{f) = ip{r), we immediately find the operator eigenvalues 



A„ = exp {2-Kin/M) , n = 0, . . . , M - 1, 



(47) 



where i here, and onwards in the similar combinations, stands for the imaginary unit i = \/— T. 

To proceed with corresponding set of constants, as a calculation trick, we assume an observation point to be 
arbitrary but close to the center of the molecule i.e. r <C i?o such that we can do approximation e~'^l'^~^^ l « e""^" 
and write the following expansions 



f- R, 



2^-^-Pi(cos7j), i = i, 



1=0 ^0 



(48) 



r ■ R 

where cos 7, = 7—; — —. 

' \r\ ■ Ro 

The addition theorem for spherical harmonics (e.g. 0) allows to write 



Pi (cos 7-,) 



47r X - 

21 + 1 



{-ly^-Yne, </>)-yr'"(7r/2, 0).e 



(49) 



Substitution of (gl]), (gT!)-!!!!]) into the eigenvalue problem 



yields 



M 

m j — 1 



M 



that is equivalent to 



C'm+1 — Ci 



This allows to conclude 
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and, consequently, 



Cj+i ^ e-'^^Ci, j = l, ...,M-1, n = 0, ...,M-1. 



Feeding this into (|46)) solves the problem of finding energy levels 



(50) 



= 0, . . . , [M/2] 



(51) 



by treating transcendental equation with respect to k„ (by means of numerical methods, e.g. the Newton-Raphson 
algorithm) 



M 
3 = i 



(52) 



Structure of the spectrum is as it follows. Ground state n — and the state with highest energy n — M/2 (in 
case of even number of atoms M) are the only non-degenerate ones. 



M 

V'o(r)=Ci^.g(|r-i?, 



, 1^ h 



3 = 1 



(53) 



M 



V,M/2(r) = Ci^(-l)^"-^g( 

3=1 



whilst all the other states are doubly-degenerate with complex conjugated wave functions 



M 



r- R, 



, K 



(54) 



(55) 



M 



^^(r) = Ci^e^(^-i)^ - .g 



f- R, 



1 i 1 



1, . . . , [M/2] - 1, 



where the constant Ci can be found from the normalization condition. 



(56) 



4.2 Continuous spectrum problem 

Now we consider the scattering of plane wave incident axially on the A/-atom potential plane symmetrical structure 
under question. 

We write the solution to the Sclirodinger equation as 



M 



V'(f) = Aoe^'=^+^C,..g 



r-R, 



(57) 



where k 



Aq is an amplitude of the incident plane wave and g (r, k) is as introduced above according to 
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As before (and using previously described notation), at each atom the solution must satisfy the ZRP condition 

m 



9 log ^ 


r- R, 




d 







(Ao + . . . + Cj_i • g {ARi, -ik) + Cj ■ ik + Cj+i ■ g {ARi, -ik) + ...) = -(3. 



\r-R,\=0 



This results in the set of linear equations 

^ /3 + ifc g{ARi,-ik) ... 
g (ARi, -ik) 13 + ik 



g{AR2,~tk) g{ARi,-ik) 
g(Ai?3, -z/c) g{AR2,-ik) 



Co 



g{AR2,~ik) g{AR3,-ik) ... (3 + ik g{ARi,-ik) 
yg{ARi,~ik) g{AR2,-ik) ... g{ARi,~ik) 13 + ik 



Cm-1 



that obviously has the solution corresponding to trivial symmetry of the system 



-^0 



1 
1 

V ) 



Ci — . . . — Cm — Com- 
From any of the equation of the set it follows that 

^0 



Com — \ 



f3 + ik + 2Y}'l;'^'^9{AR,, -iky 
Ao 



^ P + ik + 2 Ej Ji 9 (Ai?j, -^k) + g (Ai?M/2, 



M - odd, 
M — even. 



This can be unified as 



Com — — 



Ao 



l3 + ik + 2 ^(f +xm)/2-i g (Ai?^., -^k) + XM-i ■ g (Ai?A//2, ~^k) 



where xm 



M/2-[M/2] 



is an indicator of parity of the number of atoms M. 



M/2-[{M~l)/2] 

It is also useful to recall (|^ and rewrite the sum explicitly 



(M+xm)/2-1 

E 



(M+xm)/2 



exp {2ikRo sin {(j)j/2)) 
sin {4>j/2) 



with 0, = ^ (j - 1). 



Having obtained that, and taking into account the expansion 



(58) 



(59) 



(60) 



— Rj ~ r — ^ • i?j at r — >■ cx3, we can present 
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the solution (|57p as 



(61) 



where the scattering amphtude at long distances in direction n — r/r — (sin0cos(/), sin0sin(/>, cos^)"^ is given by 

M M 

0) = CoA/Mo^e-^^"-^^ =CoA/Mo5]exp{-zfci?osin0cos((/.-0j)}. (62) 
3=1 i=i 

From here the differential cross-section can be determined 



(63) 



as well as the total scattering cross-section 



f27T pTT 

^ = / / 1/(^,0)1' sin 
10 Jo 



(64) 



This also can be expressed according to the optical theorem 



a = |^Im{/(0, 0)} 



(65) 



Therefore, 



a = — — Im{CoMMo} = 



k + P 



k [k + Pf + (/3 + QY 



(66) 



where we denote 



P^^. y sm(2fci?oSm(7r,/M))^XM^^^.^^^^^^^ 



sin (TTj'/Af ) 



(67) 



1 

i?o 



(M+xm)/2-1 

E 



cos(2A:i?osin(7rj/Af)) xm- 



i • cos (2fci?o) 



sin {■nj/M) 2 
We note that for the case of a single ZRP scatterer (i.e. M = 1) characterized by 



(68) 



dlog (r'0(r)) 
dr 



r=0 



(69) 
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the formula ([5^ should be obviously replaced with 



An 



l3o + ik 



(70) 



resulting in the following formula for total scattering cross-section 



47r 47r 



(71) 



As it was announced, we believe that a simplified theoretical model of molecular structures can be built by 
replacing each multiatomic molecule with some ZRP once its spectrum is known. Indeed, in particular for scattering 
on one molecule instead of (j66p one would attempt to use (j7ip with (3q to be chosen from the discrete spectrum 
levels according to the obtained formula (HI]). 

4.3 Results for benzene molecule 

Focusing on particular case of a benzene molecule CqHq that can be considered as a cyclic structures formed by 
C — H constituents described by ZRPs, we set M = 6 and start with the bounded state problem solution. 

Structure of the discrete spectrum is as described by i.e. there are 4 energy levels with 2 middle 

ones being doubly degenerate. Numerical solution of the transcendental equations ([5^ yields the values of energy 
levels which are compared with those obtained by the well-known Huckel method on the Fig. |31 The spectrum, 
namely distances between the energy levels, reveals good correspondence of ZRP approach to the qualitative Huckel 
calculations dictating energies from the ground to the highest level to be i?o — 271, Eq — A, + A, Eq + 2A, 
respectively |17[ [5]. 



Spectra Comparison for Benzene Molecule [eV] 
ZRP calculations Huckel mettiod 





Fig. 3: Comparison of ZRP and Huckel calculations of energy levels for benzene molecule 
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To analyze continuous spectrum, we consider total scattering cross-section dependence on energy. 

We compare the results given by the standard (atomic) ZRP model of the molecule ([55)) and those followed from 
single ZRP description of the molecule, i.e. (ffTj) using the inverse of effective scattering length as introduced by 
(j4ip . As one can observe on the plots given below for differently estimated values of radius of benzene molecule 
(e.g. see the introduction of [3T]), there is moderate coincidence between conventional set of atomic ZRPs model 
and proposed single molecular ZRP approach. 



Scattering Cross Seclion lor Benzene Molecule 



Scattering Cross Section for Benzene Molecule 



Fig. 4: Total scattering cross-sections for atomic 
and molecular ZRPs: smaller radius i?oi of 
benzene molecule [C — C bond length esti- 
mation) . 



Fig. 5: Total scattering cross-sections for atomic 
and molecular ZRPs: bigger radius i?02 of 
benzene molecule (effective upper estima- 
tion). 



Also we notice that better correspondence of the molecular ZRP results to those obtained with the set of standard 
ZRPs can be achieved if in the dressing procedure one includes only the ground state and the state with maximal 
energy, i.e. Ej — — j — 0, 3. This can be viewed as if we follow the described general procedure of application 
DT for every level, but for each doubly degenerated level with Ej ~ — ^jf-, j — 1, 2 we perform an additional DT 
again with the prop function of type ()25|) but choosing a negative transformation parameter k = —Kj, j = 1,2 
(corresponding to the same energy level). This results in cancellation of consecutive contributions of the DTs with 
opposite to each other spectral parameters, hence the expression (PT|) gives 

P^{k^' + k^Y'^ (72) 

This result is not surprising since, as it can be seen from the scattering matrix formalism, such pair of DTs yields 
consecutive adding and removal of a pole on the imaginary positive k half-axis (which corresponds to a bounded 
state) of (ED): so • • 1^ = .o, J = 1, 2. 

That is effective characteristic involves only the lowest and the highest energy level parameters. In case of odd 
number of atoms this would involve only the ground state energy parameter kq. 
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Refined in this way molecular potential gives results demonstrated on the Fig. [SI [T] 




Fig. 6: Total scattering cross-sections for atomic 
and molecular ZRPs: 

smaller radius _Roi of benzene molecule, level 
degeneracy is taken into account. 



Fig. 7: Total scattering cross-sections for atomic 
and molecular ZRPs: 

bigger radius i?02 of benzene molecule, level 
degeneracy is taken into account. 



Below is the table with values of parameters used for calculations of data presented on the Fig. |3]-[71 







A 




Ro2 


1.5 Ang"^ 


-7.1 eV 


0.8 eV 


1.4 Ang 


3.6 Ang 



Table 1. Set of parameters for data plotted on the figures. 



5 Conclusions 

The presented method of extension of conventional atomic ZRP approach is made to describe molecular structures 
based on their discrete spectrum. Bounded state energies, either found from conventional ZRP calculations or in 
whatever alternative way, determine effective characteristic of newly introduced single ZRP that is expected to 
adequately describe properties of a molecule. 

Certainly, replacement of arbitrary molecules with spherically symmetric point-centers is a crude approximation 
valid only in low-energy limit, however, as it was demonstrated, even in case of highly anisotropic molecules, 
like benzene, the model gives not perfect but reasonable quality results at the same time allowing significant 
simplification in formulation of more complicated molecular structures problems with molecules instead of atoms 
as ZRP constituents. Although we believe that the molecular ZRP method should be tested on more examples 
in order to prove its quantitative merit, it seems to be a beneficial tool at least to obtain estimates of results on 
physically reasonable analytical grounds. 
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